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Abstract: Coherent QCD radiation in the hadroproduction of top quark pairs leads 
to a forward-backward asymmetry that grows more negative with increasing transverse 
momentum of the pair. This feature is present in Monte Carlo event generators with 
coherent parton showering, even though the production process is treated at leading order 
and has no intrinsic asymmetry before showering. In addition, depending on the treatment 
of recoils, showering can produce a positive contribution to the inclusive asymmetry. We 
explain the origin of these features, compare them in fixed-order calculations and the 
HERWIG++, Pythia and Sherpa event generators, and discuss their implications. 
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1. Introduction 

The observation of a substantial forward-backward asymmetry in the production of top 
quark pairs at the Tevatron [1-5] has prompted renewed theoretical study of Standard 
Model predictions for this quantity: see Refs. [6-15] and references therein. Quantities of 
particular interest are the distributions of the asymmetry with respect to some observable, 
generally denoted by O: 
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(1.1) 



where the rapidity difference is defined as Ay = yt — yi- Examples of possible observables 
O are the invariant mass m t i and transverse momentum pxtt °f t ne pair. For the former, 
QCD predicts an asymmetry that increases with m t i [6,9]. For the latter, the prediction 
changes sign as px tt increases, since the NLO loop contribution is positive at Ptu = while 
the real-emission contribution at Pxtt > is negative. After matching to parton showers 
using the MConlo prescription [16, 17], the predicted cross-over is around Pxtt ~ 25 GeV: 
see Fig. 1. 
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Figure 1: Top quark forward-backward asymmetry at the Tevatron, as a function of the transverse 
momentum of the pair: results from MC@NLO and PYTHIA. Figure from Ref. [4]. 



A surprising fact, also shown in Fig. 1, is that a leading-order parton shower event 
generator such as PYTHIA, with appropriate settings, displays a qualitatively similar p?- 
dependent asymmetry, even though the LO production processes have no asymmetry. As 
we shall see, the same is true of the HERWIG++ and Sherpa event generators, although 
their quantitative predictions differ. 1 

The explanation is a nice illustration of the QCD coherence of parton showering. It 
is true that the inclusive asymmetry built into the LO generators is zero. However, in the 
hard process qq — > ti the colour flows from the incoming quark to the top quark and from 
the antiquark to the antitop lead to a more violent acceleration of colour, and consequently 
more QCD radiation, when the top is produced backwards in the qq frame than when it goes 
forwards, as illustrated in Fig. 2. The additional radiation when the top goes backwards 
pushes the recoiling pair to higher transverse momentum. Correspondingly, events with 
forward-moving tops are left at lower transverse momentum, leading to the behaviour 
seen in Fig. 1. The effect vanishes at threshold and becomes more and more marked as 
the invariant mass of the pair increases, due to the increasing amount and scale of QCD 
radiation. 

Event generators with coherent parton showers, implemented through dipole showering 
in Sherpa and angular ordering in HERWIG++, take account of these effects. (PYTHIA uses 
a hybrid between the two.) A full NLO treatment, included in MC@NLO but not in the 
stand-alone generators, adds a finite positive virtual contribution. Nonetheless, as we 
shall demonstrate in Section 4, even LO shower models can also generate a net inclusive 
asymmetry A-p-Q > if the shower kinematics allow for migration between positive and negative 
Ay regions. 

1 For a recent review of Monte Carlo event generators, see [18]. 
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Figure 2: Colour flow and QCD radiation in (a) forward and (b) backward ti production. 

In the following section we examine in more detail the approximations made in event 
generators, in comparison to the fixed-order perturbative treatment. Then in Section 3 
we explain in general terms how they can produce a positive inclusive asymmetry while 
only containing the LO production process. In Section 4 we present results from the 
HERWIG++, Pythia and Sherpa generators for the inclusive asymmetry and various dif- 
ferential asymmetry distributions. In Section 5 we summarize our findings and comment 
on their implications. 

2. Comparison with fixed order 

To establish notation we first consider the lowest-order process, 

qipx) + <l{P2) -> 0(p 3 ) + 0(P4) , (2.1) 
for which the leading-order spin-averaged matrix element squared is 

where m is the heavy quark mass and 

s = 2pfp 2 , i = -2pi-p 3 , u = -2pt-p4. (2.3) 
The corresponding differential cross section, 



dt 167T s' 



^\M(qq^QQ)\ 2 , (2.4) 



is used for the primary hard subprocess in the event generators. Clearly, it does not 
exhibit any forward-backward asymmetry. Thus for an asymmetry to be produced by a 
leading-order generator, some parton showering must occur. 

2.1 One gluon emission 

The leading-order shower contribution is the one-gluon emission process, 

?(pi) + qipz) -> Qips) + Q(Pi) + g(k) • (2.5) 
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For the asymmetry we require the difference between this and the process 

qipx) + g(P2) -> Q(ps) + Q(P4) + • (2.6) 

The difference between the spin-averaged matrix elements squared is [19-21] 
M A = ^2\M(qq^QQg)\ 2 -^2\M(qq^QQg)\ 2 

C F (N 2 - 4) 
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(2.7) 



where 



si = (P1+P2) 2 , h = -2p 1 -p 3 , u\ = -2pi-p 4 , v\ = 2p 3 -k, 

S2 = {P3+Pi) 2 , t 2 = -2p 2 -p A , u 2 = -2p 2 -p 3 , v 2 = 2p 4 -k, (2.8) 

and Wij is the dipole radiation function 

yy.. _ Pj_\ _ / 2 .9) 

%3 \Pi-k Pj-kJ 

The asymmetries for the processes gq — > QQq and gq — > QQq are obtained from the same 
expression after crossing. They are very small and will be neglected in the following. 

We see from Eq. (2.7) that the asymmetry vanishes for N = 2. This must be the case 
to all orders, because the fundamental representation of SU(2) is pseudoreal and so in that 
case qq — > QQX is the same as qq — > QQX. 

In the event generators the gluon radiation is represented as coherent emission from the 
external lines of the Born process and so the new colour factor (iV 2 — 4)/iV is approximated 
by 2Cf = {N 2 — 1)/N. Thus we expect them to overestimate the asymmetry at non-zero 
Pt = Pt qq by around 60% in lowest order. They also neglect the second term in the square 
bracket in Eq. (2.7), which is less singular at small k, and approximate the remaining terms 
by the Born term times dipole-like factors. Thus they effectively treat the asymmetry in 
the soft gluon limit. They treat the gluon radiation more accurately in the collinear regions, 
but those regions are not dominant in the asymmetry. 

2.2 Soft gluon limit 

As explained above, we expect the event generators to reproduce the soft gluon limit of 
the asymmetry, apart from the colour factor mentioned earlier. In the limit of small k we 
have si, S2 — > s etc. and Eq. (2.7) takes the simple form 

6 (jV 2 -l)(iV 2 -4) f t 2 + u 2 2m 2 



N 3 V s 2 

x (W 13 + W 2i - W u ~ W 23 ) . (2.10) 
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In the soft gluon limit we can write the differential cross section for emission of a gluon 
with energy u as 

d 2 a oj 



J^2\ M (lQ^QQ9)\ 2 dQ , (2.11) 



dtduj (4vr) 4 s 2 

where f2 is the solid angle for soft gluon emission, and we can use the fact that [22] 
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4vr 

Fij, (2.12) 



w 2 



where fj,- is the relative velocity of i and j: 



i j 



l-(^) 2 . (2.13) 
\Pi-Pj 



Now we define the asymmetry cross section da a = da — da (3 -H- 4): 



= ^ (JV-4) ( + Fm _ ^ _ } ^ 

The radiation functions i^j appearing in Eq. (2.14) have collinear divergences along the 
beam directions, which cancel in the full expression. Regulating them with a small light 
quark mass /i, we have in the limit fi — > 0: 

F 21og(^^)-2, (2.15) 

so that 

d 2 a A 4al(N 2 -1)(N 2 -4) ft 2 + u 2 2m 2 \ ft 



dtdu, *u N 3 V s 2 + — J l0g UJ • (2 ' 16) 

In lowest order the QQ pair recoils against the emitted gluon, so the transverse mo- 
mentum of the pair is pr = — &t where dkr/kr = duo/us. We can also express Eq. (2.16) 
in terms of the qq — > QQ scattering angle 9 since 

t = -s-u = -~s(l-0cos0) , (2.17) 

where f3 = y/l — 4m 2 /s is the heavy quark cm. velocity in the Born process. Thus 

*** - "i * (Ar2 " 1){N2 ~ 4) <2- f ^ o) log ( \=£mat) . (2 . 18) 



dp T dcosd sp T N 3 V / yi + ^cos^ 

Finally, we can integrate this expression to find the recoil distribution of the asymmetry 
in the soft limit. It is convenient to normalize this to the qq — > QQ Born cross section, 2 
obtained by integrating Eq. (2.4): 

*B = ^ {N ^ni-f)£-f)- (2-19) 



2 At the Tevatron, the gluon fusion contribution to the tt production cross section is small, less than 10% 
(at LO, we find 6.7% using CTEQ6L1 PDFs [23] and scales evaluated via the top quark transverse mass). 
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Figure 3: The function F((3,0) in Eq. (2.21): the full expression (solid) and the first four terms 
of its power series expansion (dashed). 



Then we can write in general 

pr da a as {N 2 - 4) 



a B dpi 



7T 



N 



F((3,pt) , 



(2.20) 



where the coupling- and colour-stripped asymmetry function F({3,pt) is given in the soft 
limit, p T < PT ax = /3 2 v^/2, by 
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(2.21) 



Thus we expect the asymmetry in the event generators to become more negative with 
increasing top pair invariant mass y^, growing linearly with cm. velocity /3 near threshold. 
The function F(f3, 0) is shown in Fig. 3. It tends to an asymptotic value of —8 log 2 — 1 = 
—6.545 as /3 — > 1. We note that the four-term power series expansion Eq. (2.21) gives a 
good approximation over a wide range of /3. 

2.3 Beyond the soft approximation 

The coupling- and colour-stripped asymmetry function F({3,pt) remains negative away 
from the soft region but decreases in magnitude as pt increases, vanishing at the phase- 
space boundary, as shown in Fig. 4. Here, as before, (3 is defined as yl — 4m 2 /si , i.e. in 
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Figure 4: The function F(/3,px), relative to its value at pt — 0. Note that p™ ax — f3 2 \fs/2. 



terms of the overall centre-of-mass energy squared s\ = s, rather than the invariant mass 
squared of the QQ pair, S2- 

To go from this function to the p-p-dependence of the asymmetry, Afb(pt), as defined 
in Eq. (1.1), we have to include the strong coupling and colour factors and normalize to the 
differential cross section at the same pr- A full leading-order calculation of the asymmetry 
using MCFM [24] (computing ti production at NLO) yields the results shown in Fig. 5. 
For this calculation, the next-to-leading order parton distribution functions of Ref. [25] 
were used, with scale equal to the top mass {m t = 172.5 GeV), but the results for the 
asymmetry are rather insensitive to these choices. 

The real-emission asymmetry cross section in the numerator of Eq. (1.1) diverges 
as pt — > 0, due to the soft divergence discussed in the previous section. However, the 
denominator, shown on the left in Fig. 5, diverges faster, as it has initial-state collinear 
singularities that cancel in the numerator, so that ^fb(pt) is driven towards zero as pr — > 0. 

In the fully inclusive asymmetry, the divergence of the real emission contribution is 
cancelled by the singular virtual correction at pt = 0, and we have 



&b Jo dpT vr N 

where f(/3) is finite. If we limit the integration to px < qr *C m, we therefore find 

oa{pt <qr) If f m do A , f m da A 

dpr — / —j — dpx 



\Jo dp T Jg T dp T 

a s (N 2 - 4) ^ 



7T N 



F(/3,0)log(^) + G((3,q T /m) 



m J 



(2.23) 
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Figure 5: Leading-order QCD predictions for the transverse-momentum distribution of the top 
quark pair, pj, (left), and for the forward-backward asymmetry as a function of pj. (right). 



where F({3,0) is given by Eq. (2.21) and G(f3,qT/m) is regular at qT = 0. Since F((3,0) ~ 
—4 (3 is negative, a cut on px < qr adds a positive contribution to the inclusive asymmetry, 
which grows logarithmically as qT is reduced. This effect can be quite significant. For 
example, in top pair production at a pair invariant mass of 450 GeV {nit = 172 GeV, 
(3 = 0.645), a cut on px < 20 GeV gives a contribution of 3.33 as ~ 35% from the 
logarithmic term. 

In a combined NLO plus parton shower treatment such as MC@NLO, the positive 
singular virtual contribution at pt = is spread out over a finite region of pt, the so-called 
Sudakov region. This leads to a cross-over in the asymmetry from positive values at low 
Pt to negative values at higher px, as depicted in Fig. 1. In MC@NLO, the finite part of 
the virtual contribution, absent from the event generators, is also included and can affect 
the position of the cross-over. 

Monte Carlo event generators with QCD coherence will not have the correct form for 
the function G in Eq. (2.23), due to the inexact treatment of hard, non-collinear emissions, 
but they should reproduce the logarithmic term, apart from the overestimate of the colour 
factor by 60% mentioned earlier. They will also display the spreading of the positive 
asymmetry over the Sudakov region at low px- 

3. Generation of an inclusive asymmetry 

Despite the coherence effect elaborated upon in the previous section, the naive expectation 
is that the asymmetry should still sum to zero when integrated over all of phase space. 3 This 

3 Or, if using MCONLO or POWHEG [26,27], that the integrated forward-backward asymmetry should sum 
to its NLO inclusive value. 
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expectation is based on the simple fact that showers are unitary (meaning real-radiation 
corrections cancel exactly against virtual-Sudakov ones), so even though they can move 
things around in phase space, they do not generate any corrections to total cross sections. 
At the most inclusive level, this is reflected in the fact that the total integrated ti cross 
section is the same before and after showering. 

However, the asymmetry is defined in terms of two separate cross sections, one com- 
puted for Ay > and the other for Ay < 0. If the shower kinematics allow any migration 
between these two regions, then unitarity no longer guarantees complete cancellation in 
each of the regions separately, leading to the possible generation of a net inclusive asym- 
metry. Formally, we can write the cross section difference that generates the integrated 
asymmetry as 

A^+- = J da LO \ Ay>Q [A + + (l-A+)(P ++ -P + _)] 

-/^ L °Ja,<o [A_ + (1-A_)(P__-P_ + )] , (3.1) 

where the first line represents events that start (at the matrix-element level, before show- 
ering) with a positive value of Ay and the second line represents events that start with a 
negative one. The terms in parentheses represent the action of the parton shower. The 
probability for no branchings to occur is represented by the Sudakov factor, A, with sub- 
script ± reflecting that the probability to radiate can be different between an event with 
positive Ay and one with negative Ay. Indeed, as shown in the preceding section, events 
with positive Ay have less phase space for emission and so are less likely to radiate. There- 
fore, in general, we have 

A+ > A_ . (3.2) 

This, however, is not by itself enough to generate an inclusive asymmetry. The second 
terms in the square brackets in Eq. (3.1) represent those events that do experience one 
or more branchings. For these events, the final top momenta, and hence possibly their 
final rapidity difference, will depend on whether and how the top momenta are modified 
by the branchings. In the present context, we do not care about the details of how this 
occurs, merely about whether it is at all possible for an event with positive Ay at the Born 
level to migrate to negative Ay after showering, and vice versa. This is represented by 

the probabilities P^ and P (. in Eq. (3.1). If the shower model preserves the rapidity 

ordering of the tops, then 

P ++ = P— = 1 and P+- = P-+ = , (3.3) 

and so the integrated asymmetry remains zero, despite the two Sudakov factors being 
different. If, on the other hand, the shower model sometimes changes the relative rapidity 
ordering of the tops, for instance as a consequence of longitudinal recoil effects (as will 
be studied in more detail in the next section), then a total inclusive asymmetry can be 
generated. In the context of unitarity, this can be interpreted as due to the fact that 
unitarity involves an integral over the entire phase space, and hence the exact cancellation 
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that occurs in the total inclusive cross section is here broken by the splitting-up of the real- 
radiation phase space into two regions that enter with different signs in the asymmetry. 
From unitarity of the shower, we have 

P++ = 1 - P+_ and P__ = 1 - P_+ , (3.4) 

so that Eq. (3.1) can be written as 

Aa + _ = - 2 J da LO \ Ay>Q (1 - A+)P + _ + 2 j da LO \ Ay<Q (1 - A„)P_+ , (3.5) 

where we have used 

Because 1 > A + > A_, we expect the second term on the right-hand side of Eq. (3.5) to 
dominate, giving a positive inclusive asymmetry, unless there is a compensating excess of 

P-\ over P |_. However, on rather general grounds one would not expect such an excess, 

because there is less radiation when Ay > and hence a smaller probability of recoil effects 
changing the sign of Ay. Indeed we shall see below that the treatment of recoils in shower 

generators normally leads to P |_ > P^ , enhancing the positive inclusive asymmetry due 

to the unequal Sudakov factors. 

Considering Eq. (3.5) from the viewpoint of perturbation theory, we observe that the 
factors of (1 — A±) in the integrands are 0(ctg), while P±^ are O(ctg), being the con- 
ditional probabilities that gluon emission will switch the sign of Ay, given that at least 
one emission has occurred. Thus the recoil effect in showering generates an approximate 
inclusive asymmetry that starts at 0(as), like the full perturbative calculation. The fac- 
tors of (1 — A±) provide information about the virtual contribution and the probabilities 
P-j-=p specify what fraction remains after real-virtual cancellation. Since these probabilities 
depend on the strategy for treating recoils in the shower, getting the best agreement with 
the full asymmetry at C(«s) could be a good way to optimize this strategy. 

4. Comparison between parton-shower models 

In this section, we study the asymmetries produced by the following general-purpose event 
generators: 4 HerwigH — h [28] (using angular-ordered parton showers [29]), Pythia 6 [30] 
(using both its Q 2 - and p± -ordered parton-shower models [31,32], represented by tunes D6T 
and Perugia 0, respectively), Pythia 8 [33] (using pj_-ordered parton showers [34]), and 
Sherpa [35] (using p_L-ordered dipole showers [36]). Of these, HERWIGH — h and Sherpa have 
QCD coherence built in and Pythia 6 has options with varying amounts of coherence, while 
the first ISR (initial-state radiation) emission is not subjected to coherence constraints in 
this version of Pythia 8. For both Pythia 6 and Sherpa, we include some additional 
illustrations of specific shower model variations. 

A custom-made RIVET [37] analysis was used to process the events of all generators, 
ensuring uniformity of the analysis. Between 1 and 4 million events (at least) were gen- 
erated for each model. All the generators include the leading-order qq — > ti and gg — > ti 

Specifically, we use HERWIG++ 2.5.2, PYTHIA 6.426, PYTHIA 8.162, and SHERPA 1.4.0. 
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Model 


Version 


Inclusive 


m tF /GeV 




Pr,ti/GeV 






[tune] 




< 450 


> 450 


< 50 


> 50 


rlhjHWICjri — h 




1 Q 


9 7 


D.U 


O.o 


1 A 3 


Pythia 6 


6.426 [def] 


-0.1 


-0.8 


1.2 


2.5 


-42.5 


Pythia 6 


6.426 [D6T] 


-0.2 


-1.1 


1.2 


3.2 


-43.4 


Pythia 6 


6.426 [PO] 


0.8 


0.7 


1.1 


1.8 


-8.6 


Pythia 8 


8.163 [def] 


2.5 


2.4 


2.8 


2.4 


4.8 


Sherpa 


1.4.0 [def] 


5.5 


3.5 


9.2 


8.7 


-15.4 


Sherpa 


1.3.1 [def] 


6.3 


3.3 


12.1 


9.6 


-15.8 


QCD 


LO 


6.0 


4.1 


9.3 


7.0 


-11.1 



Table 1: The forward-backward asymmetry (in %) in each shower model, and to leading 

non-trivial order. is defined in Eq. (4.1). The Pythia tunes are discussed in Section 4.2.2, and 

the shower model used in Sherpa is the CSSHOWER one. The brackets after the generator version 
number denote the parameter set (tune) used, with [def] for default parameters. The fixed-order 
predictions are from MCFM [24] used to compute tt production at NLO. 



production processes, which are showered with default settings, unless otherwise specified. 
Note that we do not study the effects of matrix-element plus parton-shower matching in 
this paper. 

4.1 Inclusive asymmetry 

The inclusive asymmetry and the asymmetries with an invariant mass or transverse- 
momentum cut, 

,j(cut) I _ -(cut) I 

^(cut) = lAy>0 U iAy<0 ^ ^ 

produced by each model are given in Tab. 1. Of these, Sherpa's CSSHOWER produces the 
largest inclusive asymmetry. We interpret this as a consequence of its initial-final dipole 
kinematics [36,40,41]; part of the longitudinal momentum of the first emitted gluon has to 
come from the recoiling top quark, changing its rapidity and allowing Ay to change sign. 
Later on, we will illustrate and discuss recoil effects in somewhat more detail, in a small 
CSSHOWER case study, see Section 4.2.1 and Appendix A. 

In HERWIG++, coherence is implemented by angular-ordered parton branching rather 
than dipole showering. Parton showers associated with each incoming or outgoing hard 
parton are generated independently in angular regions defined by the colour structure of the 
hard subprocess. The showers are then combined according to a kinematic reconstruction 

J The choice of PDF set only gives small effects (< 10%) on the asymmetry, mostly via the relative fraction 
of gluon-initiated vs. quark-initiated ti production. For completeness, HERWIG++ uses the MRSTMCal 
PDF set (i.e. the LO fit from the MRST2002 family) [38], PYTHIA 6 with Perugia uses CTEQ5L [39], and 
PYTHIA 6 with D6T, PYTHIA 8, and SHERPA all use CTEQ6L1 PDFs [23]. There is also a slight dependence 
on the choice of renormalization scale, see Appendix A. 



- 11 - 



algorithm [28] that again reflects the colour structure of the subprocess. In the case of 
qq — > ti, there are two initial-final colour-connected systems, qt and qi, as illustrated in 
Fig. 2. The resulting treatment of recoils is similar to that in dipole showering of these 
systems, with a particular prescription for sharing recoil momentum within each. The 
separate recoils of the top quark and antiquark again imply that Ay may change sign. 

In contrast, the approach to initial-final dipoles in Pythia, for both p±- and Q-ordered 
showers [31,32], is to use the other incoming parton for momentum conservation, rather 
than the recoiling top. The relative rapidity ordering of the top and the antitop is normally 
preserved by this strategy, resulting in very little net asymmetry being generated. 

Overall, the asymmetries in HERWIG++ and Sherpa are comparable to the LO per- 
turbative results 6 shown in Tab. 1, suggesting that their similar recoil strategies are not 
far from optimal. 

4.2 Asymmetry as a function of top quark observables 

We now show differential spectra da / dO for four key observables and their related forward- 
backward asymmetry distributions ^4fb(0), as defined in Eq. (1.1). The observables pre- 
sented here are the azimuthal angle Ac/) between the transverse momenta of the top and 
antitop quarks, the |Ay| distribution itself, and the transverse-momentum and invariant 
mass distributions of the it pair. We take a subset of the shower versions listed in Tab. 1 
(neglecting the Pythia 6 default and Sherpa 1.3.1 versions) and compare their predictions 
with each other. 

In Fig. 6 we show the |Ay| and Acp distributions. The da/d\Ay\ predictions are very 
similar in shape; they differ in normalization because of the spread of the total inclusive 
cross sections evaluated at LO and different scales in the various event generators. The 
asymmetry rises for larger absolute rapidity differences of the top quark rapidities. The 
large \Ay\ configurations emerge more easily in scatterings with small angle to the beam (or 
the qq axis) . This produces a positive asymmetry since the associated initial-final qt and qi 
dipoles tend to emit softer gluons, with higher rate, in forward direction. The mechanism 
is also explained more fully in Section 4.2.1. All Pythia predictions increase rather mildly, 
while those of HERWIG++ and Sherpa show a steeper (approximately linear) slope, which 
is slightly flatter, but qualitatively comparable with the recent MCFM results given for 
the acceptance corrected case, see Ref. [15]. 

The Acf) variable is a typical example of an observable separating the hard-emission 
domain from the Sudakov region, here located around large Ac/) ~ tt. We depict the 
da/dA(f> spectra in Fig. 6 reflecting the different levels of hardness produced by the different 
generators with Pythia 8 giving the highest levels. Because of the strong correlation with 
the px t i observable, Acp displays a qualitatively similar behaviour of the related asymmetry 
functions. As can be seen from the plot to the bottom right in Fig. 6, the more violent 
radiation emerging from the colour dipoles spanned with backward-moving top quarks leads 
to a negative asymmetry over a wide range of angles, except for very large angles where 

6 NLO calculations of the forward-backward asymmetry in ti + 1-jet production [7, 12] give values some 
6% above the LO result in Tab. 1, e.g. in Ref. [7] a value of —3.05% is reported for > 50 GeV. 
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Figure 6: Differential cross section (in pb) and top quark forward-backward asymmetry Ape as 
a function of the modulus of the rapidity difference Ay (upper row) and the azimuthal angle A<p 
(lower row) between the top and antitop quarks. Predictions are shown as obtained from various 
Monte Carlo event generators; errors are statistical only. 



the asymmetry turns positive as a result of the Sudakov effect. While for most predictions 
the cross-over occurs at Ac/) ~ 3, two results deviate considerably from the fairly constant 
behaviour of the asymmetry (Afb ~ —0.1) for Ac/) < 2. The Pythia 6 tune D6T and 
Pythia 8 mark these two very different ends of the low-A<^> asymmetry spectrum (with 
values of —70% and +20%, respectively). For the former, soft colour coherence effects 
are overestimated whereas, for the latter, coherence effects have not been implemented, in 
particular the initial-final qt and qi dipoles are not yet treated as such. 

Fig. 7 displays the dependence of the differential cross section and asymmetry on the 
mass m t i and transverse momentum pxtt °f the tt system. As for \ Ay\, the Sudakov region 
generating positive asymmetry contributions due to soft colour coherence (A + > A_) 
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Figure 7: Differential cross section (in pb/GeV) and top quark forward-backward asymmetry Ap^ 
as a function of the mass vn t i (upper row, in GeV) and transverse momentum p T t j (lower row, in 
GeV) of the top-antitop pair. Various event generator predictions are compared with each other, 
and all errors shown are statistical only. 

applies over the entire range of the pair mass. We can use the following equation to better 
understand this behaviour and the pair mass dependence of ^4fb : 

= mf + m 2 - + 2 Et^EtI cosh Ay — 1pT,tVT,t cos ^-4> (4-2) 

where Ej, = m 2 + "p\ and A<j) is the azimuthal angle between ptj and pj> j • It is sufficient 
to focus on the cosh Ay and cos Acj) dependence of m%. The cosh Ay term is forward- 
backward symmetric and the squared mass increases with larger absolute rapidity differ- 
ences. The cos Acf) term in Eq. (4.2) may however reduce m^, but in the hard region only. 
Consequently, the cosh Ay dependence of ^4fb directly translates into a similar mass de- 
pendence. Neglecting this for a moment, we also notice that for given Ay, it is cheaper 
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to shift the masses to larger values. This is because of the enhancement of soft emissions 
(A0 ~ 7r) causing an overall plus sign in the cos A(p term. The imbalance in the soft- 
emission rate generated by soft colour coherence between the forward and backward region 
thus produces harder mass spectra in the forward region. Taken together with the |Ay| 
dependence of the asymmetry, we conclude that the cos Acj) term induces an additional 
growth of the asymmetry with increasing pair mass and a small suppression in the low 
mass region. In cases where the |Ay| dependence of the asymmetry is almost zero, the 
same mechanism may generate slightly negative asymmetries at low mass. 

Looking at the different generator results displayed in Fig. 7, the da/dm t i spectra can 
be seen to differ less than the Pxtt spectra. Again, Pythia 8 gives the hardest distribu- 
tions. As expected, we find that the tt mass dependence of the asymmetry is determined 
predominantly by the asymmetric behaviour present in |Ay|. The pattern shown on the 
upper right of Fig. 6 repeats itself here. We also observe that small negative asymmetries 
are possible for low pair masses, e.g. as shown for the Pythia 6 tune D6T. 

Finally, on the lower right of Fig. 7 one finds the results obtained by the different 
generators for the asymmetry plotted as a function of Pt.u- We already discussed the 
characteristics of Ap^ (Prtt) throughout preceding sections; thus, recalling the discussion 
of the A<p case, the Pxti results are as expected and their interpretation is straightforward. 
Exhibiting the asymmetry in terms of Px tt naturally allows for a better discrimination in 
the hard region. We thus observe that HERWIGH — h, Pythia 6 PO and Sherpa differ in their 
description of the high-p^ tail of the asymmetry, even though they sufficiently agree in the 
low-A0 region. Sherpa predicts a slope change around 50 GeV indicating a possible return 
in ApB(PT,tt) to zero for large pt, as seen in LO (Figs. 4 and 5). In contrast, HERWIGH — h 
and Pythia 6 PO maintain their trend towards more negative asymmetries. At the other 
end of the pt spectrum, the rise of the asymmetry towards lower pt is not shown by 
Pythia 6 D6T and Pythia 8, which is compatible with the findings for Acj) ~ tt in Fig. 6. 

We have studied more observables than we are able to present here, so for a more 
comprehensive comparison we refer the interested reader to the corresponding "Top quark 
(MC)" web-pages available under mcplots . cern. ch [42]. Among other things one can 
find, for example, the distributions and forward-backward asymmetries of the transverse 
momentum of the top quark, px,ti or the rapidity y t i of the tt pair; in addition many 
observables are shown separately for the high and low pair mass or transverse-momentum 
region obtained by cutting on m t j at 450 GeV or Ptu at 50 GeV, respectively. 

Summary. We can say that the HERWIGH — h and Sherpa predictions agree fairly well 
with each other, and compare quite nicely - on a qualitative level - to the -Afb (witt) an d 
Afb(|A?/|) results given in Ref. [9]. Both models incorporate soft colour coherence on a 
compatible level, which consequently may be interpreted as the source of the agreement; 
the differences lie in details such as the treatment of recoils, the shower variables and 
the form of the splitting functions used. These differences cause deviations in the high- 
Pt asymmetry spectra. The dependence on recoil effects in Sherpa is studied in more 
detail in Section 4.2.1 below. In Pythia, soft colour coherence is accounted for on a more 
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approximate level. Although the PO tune of Pythia 6 is similar to HERWIG++ in the hard- 
emission domain, it differs in the description of the Sudakov region, yielding a milder mass 
and | Ay | dependence of the asymmetry. The Pythia 6 D6T tune and Pythia 8 exhibit 
larger differences. The dependence on the shower modelling in Pythia 6 is studied in more 
detail in Section 4.2.2 below. 

4.2.1 Dependence on recoil effects: Sherpa's CSshower 

We have argued that recoil effects play an important role in producing the asymmetries 
generated in coherent parton or dipole showering. To illuminate the mechanism further, 
we have conducted a small case study based on results obtained with Sherpa's CSSHOWER. 
Some of the details have been postponed to Appendix A in order to maintain the flow of 
the main part. 

Asymmetry enhancing longitudinal recoil effects. We have identified the unequal 
Sudakov form factors in forward and backward top production as a major source of asym- 
metry. The Sudakov imbalance (A+ > A_) emerges as a result of soft colour coher- 
ence. In addition, based on Eq. (3.5) we have seen that any net migration of the type 

P |_ > P-| > 0, from the backward to the forward Ay phase space, leads to a positive 

inclusive asymmetry and an amplification of the Sudakov or coherence effect. We now want 
to trace the origin of the migration process. 

In a simple dipole picture where a gluon emission stretches (further opens the initial 
angle of) the starting initial-final qt or qt dipole, one can easily account for Ay = Ay + e 
on average. 7 Here we denote the top quark rapidity difference before the emission (or at 
the LO generation level) by Ay, and e > expresses a small positive shift. Using a simple 
generation cut on Ay, we can test this hypothesis by counting and analyzing the events 
that end up in the backward/forward region after exclusively showering tt events produced 
at LO under the constraint ±Ay > 0. 

Fig. 8 shows the main results of this migration test; more details are compiled in 
Appendix A. We have plotted the outcomes of tt production at leading, fixed order (labelled 
"before shower") and of several CSSHOWER runs (all labelled "after ...") grouped according 
to unconstrained, Ay > ("fwd") and Ay < ("bwd") LO phase-space generation: we 
distinguish between completely showered runs and runs where showering was terminated 
after just one emission. For the latter, we only show the results obtained with the restricted 
LO phase space. Focusing on Ay distributions, we make a number of observations, broadly 
confirming the physics of the simple dipole picture suggested above: 

• showering generates a positive asymmetry, which is growing with larger \Ay\ (com- 
pare "before" and "after shower" results). 

• migrations are small, happen locally but yet across the entire Ay phase space; the 
— — > + direction and, therefore, — — > + cross-overs are favoured: (1) the "bwd" 

7 In Fig. 2 (a) and (b), this corresponds to a reduction of the scattering angle, i.e. a kick of the (anti)top 
quark in the incoming (anti)quark direction, thereby mitigating the acceleration of colour charges. 
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Figure 8: The Ay distributions for different LO generation and shower modes as predicted by 
Sherpa's CSSHOWER. Dashed lines correspond to results taken from one-emission ("1 em") showers, 
while solid lines depict those after complete showering, except for the black curve, which depicts 
the fixed LO parton-level prediction in the upper subfigure. Blue [green] lines show the outcomes 
when the LO tt phase-space generation is constrained to the forward ("fwd", Ay > 0) [backward] 
region. Each lower part contains a ratio plot using the respective list's first prediction as reference. 

generated results extend into the Ay > domain filling the deficit close to the tran- 
sition region left by "fwd" generated events that now populate larger Ay, and (2) the 
migration processes in the opposite direction are largely suppressed. 

• the largest effect already originates from the first emission, the gluon emission of 
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Figure 9: The forward-backward asymmetries versus |Aj/| (upper panel) andp^ t f (lower panel) for 
the two different recoil strategies available for the CSSHOWER implementation in Sherpa. Results 
gained by the corresponding one(two)-emission(s) showers are depicted as well. Each subfigure is 
supplemented by a ratio plot using the default CSSHOWER prediction for reference. 



the initial-final qt and qi dipoles (compare dashed with solid lines). 8 This does not 
necessarily mean that the corrections from multiple emissions are negligible; they 
easily give 10-20% effects, as can be seen in both Figs. 8 and 9. 

This pattern carries over to the high-p^ region, pj> t i > 50 GeV, except for the fact that the 



s This is also reflected by finding comparable Afb when running the CSSHOWER in one-emission mode 
(4.1%, 2.7%, 6.6%, 6.9% and -17.1%, to be compared to the SHERPA 1.4.0 entries of Tab. 1, in this order). 
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total asymmetry now turns negative, cf. the lower panel of Fig. 8. The migration is more 
severe, but cannot overcome the overall negative trend, caused by the more violently radi- 
ating "bwd" generated initial dipole configurations: the radiation imbalance predominates 
over the migration effect. 

Comparison of recoil strategies. Sherpa's CSshower implementation provides two 
recoil schemes, the default one, see Refs. [40,43] and the original CS scheme as advocated 
in Refs. [36,44,45]. They differ mainly in their treatment of the transverse recoils. The 
distribution of the longitudinal recoil momenta effectively is the same in both schemes, as 
documented in the top panel of Fig. 9 where we show Ape as a function of | Ay| (exhibiting 
all of the characteristics described earlier). 9 We see again, the bulk of the asymmetry is 
already produced by the one(two)-emission(s) showers, similarly for j4fb (Pr,tt) displayed 
in the lower panel of Fig. 9. The pT-dependent asymmetry function is the prototype of dis- 
tributions discriminating clearly between the two recoil strategies: the original CS scheme 
is more like that in the Catani-Seymour NLO calculational scheme [44]: the (transverse) 
recoil from a gluon emitted off a qt configuration (and likewise off a qi one) is compensated 
by the top quark, regardless of its role in the emission process (emitter or spectator). The 
prediction given by the original scheme therefore remains flat at about —15% for high pxtt 
while the default prediction levels off close to zero from below. In the default scheme, the 
gluon recoil is rather divided over the entire set of final-state partons. This requires an 
additional transverse boost combined with a rotation that in turn washes out the radiation 
imbalance between the forward and backward regions for very large tt- 

4.2.2 Dependence on shower model: Pythia 

For processes like ti production, which do not contain any QCD jets at the Born level, both 
Pythia 6 and Pythia 8 use so-called "power showers" [46] to populate the tij phase space. 
Since the LL splitting kernels generally represent an overestimate in the region of very hard 
jets, a factor that suppresses such emissions has been introduced in the pj_-ordered showers 
in both Pythia 6 and Pythia 8, similar in spirit to a matrix-element correction [47] but 
with a much simpler analytical structure. In Pythia 8, the suppression factor is derived 
from universal t-channel arguments [34] and does not depend on the colour structure of the 
event, wherefore it does not contribute to the generation of any ti asymmetry. In Pythia 6, 
the suppression factor is [48] 



where Pqj corresponds to the parameter PARP(67) in the code, prevoi is the evolution scale 
for the branching, and su is the invariant mass squared of the radiating parton with its 
colour partner, with all momenta crossed into the final state (i.e. it is s for annihilation- 
type colour flows and —t for an initial-final connection). This is motivated partly from 

9 The Afb (witt) are also broadly unaffected by the scheme change; similarly the ^-pg^ deviate only 
marginally from the default SHERPA 1.4.0 numbers stated in Tab. 1: 5.4%, 3.1%, 9.7%, 8.3% and -17.8%, 
to be compared in this order. 




(4.3) 
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Figure 10: Transverse-momentum distribution da/dpT } ti ( m pb/GeV) of the top quark pair (left) 
and the related top quark forward-backward asymmetry (right) in five different tunes of Pythia 6. 



studies of similar factors in the context of "smooth ordering" , introduced in [49] . Finally, 
in the Q 2 -ordered shower model in Pythia 6 [30], a veto on the emission angle is placed, 
which depends explicitly on the direction of the colour partner. 

To illustrate the effect of these choices on the asymmetry, we show the dependence of 
the asymmetry on pxtt f° r fi ve Pythia 6 tunes in Fig. 10. 

The (5 2 -ordered D6T tune includes the explicit angular cut on the emission angle 
mentioned above, which a priori should produce an effect qualitatively similar to that of 
the angular-ordered showers in HERWIGH — h However, in Pythia's Q 2 -ordered shower, the 
effect is amplified, as follows. When an initially massless ISR parton evolves to become a 
massive jet, its virtuality is generated by reducing its momentum while keeping its energy 
unchanged. The momentum of whatever the ISR parton is recoiling against (here the ti 
system) is then also reduced to ensure momentum conservation. Thus, to get a ti pair 
with a certain pr, the Q 2 -ordered shower must first radiate a massless ISR parton with 
an initially much larger value of pt, which is then reduced when that parton acquires a 
virtuality. The fact that this compresses the pt spectrum is well-known from the Drell-Yan 
case (cf. e.g. [42]) and can be counteracted by choosing a low renormalization scale for as 
(as has been done in D6T, which uses //r = 0.45 Pt), thus bringing the px spectrum itself 
back into rough agreement with other models, as illustrated in the left-hand panel of Fig. 10. 
The asymmetry spectrum, however, remains compressed, as shown in the right-hand panel. 

Among the four p^-ordered tunes, the Perugia (P0), Perugia HARD (PHARD), and 
Perugia SOFT (PSOFT) ones [48] use the suppression factor defined by Eq. (4.3), while 
Zl does not apply any suppression. The central Perugia (P0) tune uses Pqi = 1. This 
generates an asymmetry that begins to turn on at roughly p? = mt/2. The HARD and 
SOFT variations use Pqj = 4 and Pqj = 1/4, respectively, which modifies the turn-on 
point. Neither the HARD nor the Zl tunes exhibit any significant asymmetries in the 
region plotted here, similarly to the case in Pythia 8. 
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5. Summary and implications 

The studies presented above arose from the initially surprising observation that Monte 
Carlo event generators can produce non-zero forward-backward asymmetries in top pair 
production, even when treating the relevant subprocess qq — > ti at leading order, which 
has no such asymmetry. Our studies show that these asymmetries arise from valid physics 
built into generators with coherent parton or dipole showering. While not quantitatively 
correct in every detail, the coherent showering approximation captures essential features of 
the physics not hitherto well understood, which may serve as a guide to the contributions 
of higher orders. 

The generated asymmetries are of two kinds. First, in the differential cross section 
at non-zero transverse momentum of the top pair, a negative asymmetry results from the 
extra QCD radiation emitted when the top quark is produced backwards in the rest frame 
of the pair. This effect is manifest in the qq — > tig matrix element and is present in the 
generators in the soft approximation, with a colour coefficient that is exact in the large- N 
limit but 60% too large at N = 3. 

In fixed-order perturbation theory, the asymmetry at non-zero px of the pair tends to 
zero from below as pt — > 0. However, precisely at pt = there are singular virtual con- 
tributions that lead to a positive overall inclusive asymmetry which grows with increasing 
invariant mass of the pair. The event generators perform an approximate all-order re- 
summation of perturbation theory, which smears out the singular contributions at pr = 
over a finite Sudakov region, and so the asymmetry changes sign at some point and be- 
comes positive at small pr- The precise switching point is sensitive to finite terms and 
higher-order corrections not included in the generators, but the change of sign is a striking 
general prediction that should be investigated experimentally. 10 It is also worth noting 
that, because of this switch, a bias towards low pr, or against extra jet production, in 
the method used to reconstruct the tops could lead to a significant upward shift in the 
measured inclusive asymmetry. 

The other type of generated asymmetry is an overall inclusive one, positive and growing 
in value with increasing invariant mass of the pair. In fixed-order perturbation theory, such 
an asymmetry appears at order as relative to the Born process and is due to a positive 
asymmetry in the virtual correction, which dominates over the negative contribution of real 
emission discussed above. The event generators implicitly contain virtual corrections, in the 
form of the Sudakov factors that drive the showers and produce pt smearing. These factors 
are a reflection of unitarity, which implies that showering cannot change the inclusive cross 
section from the value established by the primary subprocess. One might therefore think 
that the inclusive asymmetry also could not change from zero when the primary process is 
symmetric. However, that is not the case, because the asymmetry is not an inclusive cross 
section and so is not protected by unitarity. 

In fact, the same effect that generates a negative asymmetry at non-zero pr, namely 
the extra radiation in backward top production, tends to produce a positive inclusive 

For a related quantity, the asymmetry when there is an extra jet with pr > 20 GeV, a negative value 
has indeed been reported [4]. 
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asymmetry. As expressed by Eq. (3.5), it arises from the difference between the Sudakov 
factors for forward and backward top production and the migration of recoiling top quarks 
between hemispheres. Thus it is a combination of real and virtual effects, which is of relative 
order as, because the difference of Sudakov factors is of that order while the forward and 
backward migration probabilities are pure numbers, modulo higher-order corrections, with 
magnitudes that depend on how the event generator treats recoils. The fact that recoil 
strategies based on colour flow produce inclusive asymmetries, similar to the full fixed-order 
one, suggests that the asymmetry can be regarded as arising in this way. Such a viewpoint 
could serve as a guide towards the more correct treatment of recoils, and conversely as an 
indication of the possible effects of higher orders beyond the range of explicit calculations. 11 

We believe that these findings have important implications for the interpretation of the 
experimental data. At the very least, one needs to be aware that the available event gen- 
erators can produce significant asymmetries where none were previously expected. Monte 
Carlo estimates of corrections to asymmetries could be affected by this, particularly cor- 
rections to "parton level". Moreover, these corrections will likely be model-dependent as 
documented by the detailed parton-shower comparison presented here. The results de- 
pend on the way colour coherence is implemented in the various codes, which we have 
summarized at the end of Section 4.2. 

On a more theoretical level, the fact that these asymmetries are due to recoils points 
to the importance of recoil effects, which are often neglected in estimates of higher orders 
based on soft gluon resummation. 

There are clearly many directions in which the studies presented here could be ex- 
tended. The effects of asymmetries produced or enhanced by parton showering in gen- 
erators that match to an NLO calculation, such as MCtsNLO and POWHEG, need to be 
assessed. Similarly for schemes that match to LO multi-parton matrix elements. Analo- 
gous effects will also be present in the ti charge asymmetry, currently being investigated 
at the LHC, and in other processes where the colour flow of the primary process affects 
parton showering. 
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Progress has been reported in the literature, see for example [7,12,50]. 
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Figure 11: (Left) Sherpa CSSHOWER predictions for the p T ^t distribution using different LO 
generation and shower modes. Dashed fines correspond to results taken from one-emission ("1 em") 
showers, while solid lines depict those after complete showering. Blue and green lines show the 
outcomes under constrained LO phase-space generation, to the forward ("fwd", Ay > 0) and 
backward region, respectiveiy. (Right) Sherpa CSSHOWER results for Ape (PT,ti) obtained under 
scale choice variations (mj^t-like, i.e. default versus m t ) combined with simultaneous variation of 
yitR and fj,F, which was also applied to the parton showering. 



A. Appendix: Additional Sherpa CSshower studies 

Migration tests. In the left part of Fig. 11 the radiation imbalance between "fwd" and 
"bwd" initial dipole configurations, seen in Fig. 8, is clearly documented by the respective 
predictions for the pt distribution of the tt pair. The LO configurations emerging from the 
"fwd" (Ay > 0) phase space generate a more steeply falling px spectrum with respect to 
that produced by the "bwd" (Ay < 0) generated dipoles. The turn-over already appears 
a t Pt,u ~ 10 GeV, slightly below the turn-over found in j4fb (Pt it) ? cf- Figs. 9 (bottom) 
and 11 (right). The upward shift is a result of the migration in Ay. 

Another way to express the results of Fig. 8 uses the after-showering, forward (Ay > 0) 
cross section fractions, which we define as r^^' = cr( cut )J Ay>0 / a^ cut \ Tab. 2 shows them 
for different before-showering, tt rapidity-difference regions, listed in ascending order of 
rf W( j. This is to summarize, on a more quantitative level, all earlier findings: considerably 

larger migration in > + than opposite direction (rows 1, 2 versus 4, 5); a factor ~ 3 

increased migration in both directions for harder emissions (pj< t j > 50 GeV) ; milder migra- 
tion in the low-pr with respect to the high-p^ region. The increasing pf-veto efficiencies 
reflect once more the radiation imbalance between forward and backward phase spaces. 12 
Focusing on the near-transition regions (rows 2 and 4), we observe an enhanced migra- 
tion activity with respect to that found for both hemispheres (rows 1 and 5). This nicely 
confirms the locality assumption, Ay = Ay + e. 

12 Again, the entire pattern already occurs after one emission: for the Ay < mode, the f^*' drop by 1% 
only. On similar note, the numbers barely change when switching to the other CSSHOWER recoil scheme. 
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Pr,tt/GeV > 50 



phase space 



rfwd [%] 
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fwd 



< o 

[-0.4,0] 
[-co, oo] 
[0,0.4] 
> 



5.9 
14.7 
52.7 
98.2 
99.2 



83.6 
87.8 
88.0 
90.6 
92.4 



3.7 
10.2 
54.1 
98.5 
99.3 



16.7 
46.8 
42.7 
95.2 
97.7 



Table 2: Forward (Ay > 0) cross section percentages, r^ d , according to Sherpa's CSSHOWER, 
where = (J ' cut ' ) _| Aa>0 / cr ' cut '- The LO cross section a is 4.94 pb, dropping to 0.955 pb under 

the influence of the rapidity-difference generation cuts < ±Ay < 0.4 applied to LO ti hadropro- 
duction. For the \ow-pr region, PT,tt < 50 GeV, the cut efficiencies are stated explicitly. 

Scale variations. We have checked the scale dependence of the SttERPA CSSHOWER pre- 
dictions, which we illustrate for A-pB{PTtl) m the right panel of Fig. 11. The impact of 
using the default, mj^-like scale choice (m\ t = nrf + Pj>t) versus a fixed nit scale, and 
varying the renormalization and factorization scales simultaneously in each case by factors 
of 2, is seen to be of the order of ±20% at intermediate Pxtt- This order of magnitude is 
consistent with that expected for a leading-order quantity, from variation of the scales and 
PDFs in the matrix element and parton showers, and is generally less than the variation 
due to different ways of treating recoils, or assessing the Pxtt distributions themselves. 
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